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ABSTRACT 


Two  aerodynamic  methods  of  creating  transmission 
"windows”  in  typical  plasma  sheaths  surrounding  re-entry 

vehicles  have  been  investigated:  (l)the  use  of  jets  of  coolant 
from  flush  nozzles  which  penetrate  across  the  plasma  slab 
and  then  mix  with  the  ionized  ga  sc  s  a  s  they  a  re  carried  down¬ 
stream,  (2)  the  use  of  a  physical  flow  diverter  to  bypass  some 
of  the  ionized  flow  laterally  around  the  antenna.  Equations 
are  developed  for  the  lateral  jet  penetration  in  a  supersonic 
stream  and  for  one  -dimensional  mixing  for  both  perfect  gases 
and  real  gases  in  thermodynamic  equilibrium.  The  assump¬ 
tion  of  real-air  equilibrium  thermodynamics  yields  coolant 
requirements  considerably  below  the  perfect  gas  values  for 
jet  injection  method.  It  is  concluded  that  the  jet  injection 
method  v/i  11  reduce  a  200  me.  signal  attenuation  to  3  db.  with 
n&edest  jet  flowrates  of  helium  on  a  typical  re-entry  trajec- 
t©*?y.  Conversely,  the  use  of  a  physical  flow  diverter  is  not 
an  effective  method  of  creating  a  transmission  window  where 
tbs^ionization layer  is  not  adjacent  to  the  vehicle  surface  or 

where  the  attenuation  needs  to  be  reduced  more  than  10  db. 
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FOREWORD 


This  report  was  prepared  by  MITHRAS,  Inc., 
of  Cambridge,  Massachusetts,  for  the  Microwave 
Physics  Laboratory,  Air  Force  Cambridge  Research 
Laboratories,  Hanscom  Field,  Bedford,  Massachusetts, 
under  Contract  AF  19(628)- 3852.  The  work  was  initiated 
by  Mr.  Walter  Rothman  and  monitored  by  him  and  Mr. 
John  F.  Lennon. 

The  investigations  whose  results  are  reported 
here  were  conducted  during  the  period  14  October  1963 
to  14  October  1964.  They  were  directed  by  Mr.  Jacques 
A.  F.  Hill  and  carried  out  by  Mr.-$q£eph  J.  Jtossi  and 
Mr.  Richard  H.  Adams.  This  report  was  written  by 
Mr.  Joseph  J.  Rossi  and  Mr.  Richard  H.  Adams.  The 
research  program  is  being  continued  under  the  same 
contract  number. 
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1.  INTRODUCTION 


The  plasma  sheath  surrounding  a  hypersonic  re-entry  vehicle 
attenuates  the  radio  signals  transmitted  to  and  from  the  vehicle,  often 
to  the  point  of  complete  blackout.  Much  work  has  been  done  on  the 
physics  of  the  interaction  of  the  plasma  and  the  electromagnetic  wave 
near  the  vehicle,  and  experimental  flight  data  supporting  theories  and 
prediction  is  now  becoming  available. 

On  experimental  manned  re-entry  vehicles  the  blackout  problem 
has  been  solved  by  on  board  recorders  which  store  data  during  the  black¬ 
out  and  then  playback  the  data  once  the  blackout  region  has  been  passed. 
The  manned  Mercury  vehicles  have  merely  “lived  with11  the  blackout 
problem.  However,  with  the  advent  of  the  Apollo  program  and  other 
manned  vehicle  programs, particularly  manned  maneuvering  vehicles, 
the  importance  of  uninterrupted  radio  contact  with  the  vehicle  during  re¬ 
entry  has  increased.  This  necessitates  developing  methods  of  creating 
“windows"  through  the  plasma  sheath. 

Under  contract  AF  19  (628) -3852  with  the  Microwave  Fhysics 
Laboratory  of  AFCRL,  MITHRAS,  Inc.  ,  has  conducted  a  study  of  aero=- 
dynamic  methods  of  creating  radio  windows  through  the  plasma  sheath 
surroun-Hng  a  re-entry  vehicle.  In  the  course  of  this  study,  two  techniques 
have  been  developed.  One  uses  jets  of  coolant  from  flush  nozzles  which 
penetrate  across  the  layer  of  ionized  flow  and  then  mix  with  it  as  they  are 
carried  downstream.  The  other  uses  a  physical  flow  diverter  to  bypass 
some  of  the  ionized  flow  laterally  around  the  antenna.  The  region  of 
separated  flow  behind  the  flow  diverter,  through  which  the  antenna  looks, 
is  kept  cool  by  low-velocity  coolant  injection. 
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2.  SUMMARY  AND  CONCLUSIONS 


The  plasma  sheath  properties  on  typical  re-entry  vehicle  shapes 
have  been  investigated  by  treating  the  sheath  as  a  homogeneous  slab 
bounded  by  free  space.  Typical  shapes  for  which  the  attenuation  of  a 
200  Me.  signal  was  determined  during  re-entry  are  the  cone,  the  flat 
plate,  and  a  blunt  nosed  body. 

Two  aerodynamic  methods  of  creating  transmission  '’windows"  in 
the  plasma  sheath  have  been  investigated;  (1)  the  use  of  jets  of  coolant 
from  flush  nozzles  which  penetrate  across  the  plasma  slab  and  then  mix 
with  it  as  they  are  carried  downstream,  (2)  the  use  of  a  physical  flow 
diverter  to  bypass  some  of  the  ionized  flow  laterally  around  the  antenna. 
Equations  were  developed  for  the  lateral  jet  penetration  in  a  supersonic 
stream  and  for  one-dimensional  mixing  for  both  perfect  gases  and  real 
gases  in  thermodynamic  equilibrium.  Charts  of  coolant  required  to  achieve 
transmission  of  a  200  Me.  signal  on  typical  vehicles  were  derived. 

Finally,  the  effect  of  the  use  of  a  side  jet  on  the  vehicle  aerodynam¬ 
ics  was  briefly  investigated. 

The  results  of  the  above  summarized  investigations  have  led  to  the 
following  conclusions  with  regard  to  the  problem  of  achieving  a  transmission 
"window"  in  a  plasma  sheath. 

(1)  With  the  approximation  of  thermodynamic  equilibrium, 
transmission  losses  at  200  Me  ore  negligible  for  tem¬ 
peratures  below  3000PK  at  all  pressures  of  interest. 

(2)  If  it  is  assumed  that  an  average  transmitting  frequency 
for  re-entry  vehicles  is  in  the  200  megacycles  range 
and  the  plasma  thickness  is  0.2  ft.  ,  then  transmission 
is  possible, 

(a)  for  a  cone  only  if  M  sin  0  <6 
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(b)  for  a  flat  plate  only  if  the  angle 

of  attack  and  wing  loadings  are  low 

(c)  for  a  blunt  body  only  far  back  on 
the  body. 

(3)  Jet  penetration  across  the  entire  shock  layer  at 
typical  antenna  installations  maybe  achieved  with 
practical  nozzle  geometries  and  coolant  supply 
pressures . 

(4)  The  jet  method  can  reduce  the  plasma  temperature 
to  levels  at  which  transmission  is  possible  with 
modest  jet  mass  flows. 

(5)  The  assumption  of  real-air  equilibrium  thermo¬ 
dynamics  yields  coolant  requirements  considerably 
below  the  perfect-gas  values  for  the  jet  injection 
method. 

(6}  The  effects  of  the  side  jet  on  the  re-entry  vehicle 
aerodynamics  are  slight,  typically  causing  a  trim 
angle -of-attack  less  than  one  degree. 

(7)  The  use  of  a  physicalflow  diverter  is  not  an  effective 
method  of  creating  a  transmission  window  where  the 
ionization  layer  is  not  adjacent  to  the  vehicle  surface 
or  where  the  attenuation  needs  to  be  reduced  more 
than  10  db. 
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3.  PLASMA  SHEATH  PROPERTIES 
ON  TYPICAL  VEHICLES 


3.  1  Transmission  Through  a  Plasma  Sheath 


The  electromagnetic  wave  incident  on  a  boundary  such  as  the 
discontinuity  in  going  from  an  antenna  surface  to  an  ionized  flow  field 
reflects  some  energy  at  the  surface  and  transmits  some  energy.  This 
problem  may  be  simplified  by  considering  only  normally  incident  waves. 

If  the  plasma  slab  is  considered  bounded  on  both  sides  by  free  space  and 
considered  homogenous  rather  than  containing  the  bightemperature  and 
density  gradients  generally  surrounding  a  hypervelocity  vehicle,  further 
mathematical  simplifications  are  possible. 

When  an  electromagnetic  wave  is  imposed  upon  a  medium  con¬ 
taining  free  electrons  such  as  a  plasma,  each  electron  may  be  considered 
to  act  as  a  miniature  antenna  and  to  absorb  energy  from  the  wave.  Since 
the  electron  is  accelerated  by  the  electric  field,  the  energy  of  the  wave 
may  be  thought  of  as  being  converted  into  energy  of  the  electrons.  This 
conversion  of  energy  results  in  a  decrease  of  signal  strength  through  re¬ 
flection  and  absorption  effects.  A  tractable  expression  for  the  attenuation 
in  transmitting  through  a  plasma  sheath  can  be  derived  by  solving  Maxsvell1  s 
equations  and  the  equations  of  motion  for  the  electrons  in  this  field.  Plane- 
wave  solutions  to  Maxwell’s  equations  for  homogenous  media,  for  the  case 
of  a  bounded  plasma  shown  in  Figure  3.1,  were  presented  in  Reference  1. 

This  theory  was  used  to  calculate  transmission  losses  through  a  plasma 
sheath  since  analyses  involving  more  mathematical  rigor  are  extremely 
laborious  and  many  times  present  insolvable  expressions.  In  view  of  other 
related  uncertainties  (i.e.  non- equilibrium  effects,  viscous  interactions, 
and  electron-heavy  particle  collision  cross- section)  such  rigor  is  not  clearly 
warranted.  For  example  Figure  3.2,  taken  from  Reference  2,  shows  approxi¬ 
mately  a  35  db  difference  in  transmitting  through  a  bounded  plasma  slab 
when  a  non-equilibrium  rather  than  an  equilibrium  gas  flow  field  is  assumed. 


4 


MC  63-78-R1 


This  change  due  to  the  selected  gas  model  is  significantly  greater  than 
what  one  would  expect  if  more  complex  plasma  models  were  chosen. 

The  solutions  for  a  bounded  plasma  slab  in  terms  of  signal 
attenuation  as  a  function  of  plasma  frequency,  ,  collision  frequency, 
and  signal  frequency,  to,  for  thickness,  6,  of  0.  1  ft.  and  0.2  ft,  are 
presented  in  Figure  3.3.  This  figure  illustrates  the  relative  impor¬ 
tance  of  the  various  parameters  over  the  range  of  interest  for  re-entry 
vehicles.  The  3  db  attenuation  line  was  chosen  as  a  reasonable  blackout 
limit, 

3.  2  Plasma  Sheath  Properties  Related  to  Flow 

The  solutions  shown  in  Figure  3.3  for  a  plane  vrave  in  passing 
through  a  plasma  slab,  if  an  equilibrium  gas  is  assumed,  show  that  the 
attenuation  of  a  signal  is  afunction  onlyofthe  following  four  variables; 

to  signal  frequency 

Up  plasma  frequency 

toc  electron  collision  frequency 

6  plasma  thickness 

In  turn,  in  the  case  of  an  equilibrium  gas,  both  the  plasma  frequency  and 
the  collision  frequency  are  solely  a  function  of  the  state  properties  of  the 
gas.  Once  these  properties  are  known  the  attenuation  resulting  in  trans¬ 
mitting  through  a  plasma  sheath  with  a  given  signal  frequency  and  plasma 
thickness  can  be  solved  for. 
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The  relations  describing  the  olasma  and  collision  freauencies  are 

-  - - - p  '  *-  —  -  * 

as  folllows: 


to 

c 


c 


N. 

J 

a 

j 


56.4  JnH*  RAD/SEC 
3 

electron/  meter 

c  Z  N.  Q. 

J  J 

average  velocity  of  electron 
concentration  of  j**1  species 

electron  cross-section  of  j**1  species 


(3.1) 


(3.2) 


The  estimate  for  the  electron  cross-section  has  been  taken  from  Refer¬ 
ence  3  and  results  in  the  following  formula  for  the  collision  frequency: 

to  =  5.0  •  1012  p  (3.3) 

c  x 

T  =  °  K 

P  =  ATM 


The  electron  density  for  real  air  can  be  determined  from  Feldman 

(4) 

Lin's  curves  ,  however,  can  be  approximated  to  within  a  factor  of  2  by 
the  following  empirical  formula  from  Reference  3. 


N 

e 


1. 87  •  10 


20 


0.6 


( 


T  -  1400 
3500 


8.  7 


(3.4) 


1800  <  T°K  <  7000 


10"6  <  ^P/P0^  <  10 
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JL  VI 


4-1 - 


vaic 


log 


10 


T  °  K 
6000 


0.6  (3.5) 


By  using  equations  3.  1  through  3.  5  plus  the  relations  describing  the 
attenuation  of  a  signal  passing  through  a  bounded  plasma  slab  with  the 
assumption  of  a  constant  signal  frequency  and  plasma  thickness,  the 
attenuation  may  be  found  as  a  function  of  pressure  and  temperature. 

For  illustrative  purpose  a  signal  frequency  of  200  megacycles  and 
plasma  thickness  of  0.  1  ft.  and  0.  2  ft.  were  chosen  for  the  results  pre¬ 
sented  in  Figure  3.4  through  3.  7,  and  throughout  this  study.  These  figures 
describe  the  temperature  criteria  necessary  for  cooling  sufficiently  to  pre¬ 
vent  blackout.  It  is  seen  that,  if  the  local  flow  over  an  antenna  is  cooled  to 
approximately  3000°  K,  and  equilibrium  is  achieved,  transmission  is  clearly 
possible.  The  use  of  electrophilic  gases  in  combination  with  an  injected 
coolant  would  allow  a  further  decrease  in  atten\_ation  which  is  not  shown  in 
these  figures. 

3.  3  Example  of  Basic  Configurations 

In  order  to  establish  the  magnitude  of  the  problems  of  transmitting 
through  a  plasma  sheath,  several  basic  configuration  shapes  were  chosen 
for  comparison.  No  specific  vehicle  has  been  chosen  for  analysis  since 
it  is  felt  that  a  flexibility  of  choice  of  the  vehicle  should  be  maintained  so  that 
the  results  of  the  study  may  have  the  broadest  application  possible.  With 
this  in  mind,  a  few  basic  vehicle  shapes  which  could  in  composition  make 
up  a  more  complex  vehicle  shape  were  studied.  The  shapes  studied  were 
a  cone,  a  flat  plate  at  angle -of- attack  and  a  blunt  body. 

The  totai  attenuation  of  a  200  Me.  signal  transmitted  through  a  0.  1 
ft.  and  0.  2  ft.  plasma  slab  on  these  various  bodies  was  taken  as  a  realistic 
example  of  interest  in  the  study.  As  will  become  evident  in  section  4.3,  the 
local  velocity  in  the  region  of  the  plasma  is  a  most  significant  parameter  in 


7 


MC  63-73-R1 


the  injection  of  a  coolant  to  reduce  attenuation. 


3.  3.  1  Cones 


For  the  data  dealing  with  cones,  no  specific  trajectories  are 
required  as  the  plasma  frequency  on  cones  may  conveniently  be  related  to 
the  hypersonic  similarity  parameter,  M^'sin  0c>  as  suggested  by  Romig 
(Reference  5). 

The  various  flow  properties  on  the  surface  of  a  cone  as  related  to 
sin  ©c  are  presented  in  Figures  3.  8  through  3.  11.  The  density  and 
temperature  plots  (Figures  3.  8  and  3.  9)  along  with  the  relations  of  section 
3.2  allow  the  calculation  of  plasma  frequency  and  collision  frequency  as  pre¬ 
sented  in  Figure  3.12.  The  shock  layer  thickness  may  be  determined  from 
the  shock  wave  angle  as  defined  inFigure  3.10.  The  cone  surface  velocity 
presented  in  Figure  3.11,  as  pointed  out  previously,  is  an  important  para¬ 
meter  when  the  injection  of  a  coolant  is  considered. 

The  relationship  of  total  attenuation  of  a  200  Me  signal  through  a  0.1 
ft  and  0.2  ft  bounded  plasma  for  a  cone  at  an  altitude  of  210,000  ft  as  a 
function  of  the  hypersonic  similarity  parameter  is  presented  in  Figure  3.13. 
This  figure  was  derived  from  the  plasma  and  collision  frequency  data  of 
Figure  3.  12,  and  used  in  conjunction  with  the  bounded  plasma  attenuation 
data  of  Figure  3.3.  Note  the  close  relationship  between  the  shape  of  the  total 
attenuation  curve  of  Figure  3.13  and  the  shape  of  the  plasma  frequency  curve 
of  Figure  3.12.  Bearingthis  similarity  in  mind,  and  observing  that  transmis¬ 
sion  is  possible  at  210,  000  ft  for  cones  if  sin  8^  <  10,  then  it  may  be  con¬ 
cluded  that  transmission  of  a  200  Me  signal  is  possible  down  to  55,  000  ft 

altitude  for  cones  of  plasma  thickness  on  the  order  of  0.  2  ft  and  M  sin  8  ^6 

oc  c 

3.3.2  Flat  Plate 


Data  computed  for  a  flat  plate  was  derived  by  assuming  that  a 
body,  one  of  whose  surfaces  was  a  flat  plate,  followed  an  equilibrium  glide 
trajectory  at  an  angle  of  attack  for  maximum  C^.  This  configuration  was 
chosen  as  the  two-dimensional  flat  plate  pressures  are  applicable  to  the  under 
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surface  of  a  flat-bottomed  body,  delta  wing,  etc.  For  computational 
purposes  the  shock  wave  angle  was  assumed  constant  at  53°  44*.  The 
calculated  trajectories  for  different  values  of  W/C^A  are  shown  in 
Figure  3.14  and  the  surface  velocities  on  the  plate  are  presented  in 
Figure  3.15. 

*",e  plasma  frequencies  for  varying  W/C^A  are  presented  in 
Figure  3.16.  Figure  3.17  shows  the  trend  of  plasma  and  collision  fre¬ 
quencies  with  W/C^A  for  the  specific  altitude  of  210,000  ft.  The  attenua^- 
tion  as  a  function  of  W/C  A  for  210,000  ft  altitude,  a  transmission  fre- 
quencyof200  Me  and  plasma  thickness  of  0.1  ft  and  0.2  ft  is  presented 
in  Figure  3.18.  Note  the  similarity  between  this  figure  and  Figure  3.13 
which  presented  the  attenuation  on  a  cone  for  the  same  conditions.  The 
strong  function  of  attenuation  with  W/C^Ais  evident  in  Figure  3.18.  A 
low  W/CT  A,  or  a  lowwing  loading  lowers  the  attenuation  rapidly  on  a  flat 

l~i 

plate  as  would  be  expected  since  this  has  the  effect  of  dissipating  the  kinetic 
energy  of  the  body  in  a  region  of  low  density  and  therefore  small  electron 
concentration. 

3.  3.  3  Blunt  Body 

Pressures  on  a  blunted,  hemisphere  cylinder  are  given 
by  the  following  equation  from  Reference  6: 


P 

x 

P~ 

00 


£  (V)  M«2  CD 
x/d 


i/2 


+  0.44 


This  equation  has  been  derived  for  equilibrium- air  at,  speeds  consistent 
with  b?  *st  wave  theory.  The  formula  has  been  compared  with  experimental 
equilibrium  gas  flows  and  found  to  be  in  good  agreement  (Reference  6). 

Px  is  surface  pressure  along  body 
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•ivy; 


is  cl  tabulated  function  Ox  the  iSeuiropic 

exponent  at  the  stagnation  region 


S 

X/D 


=  0.097  drag  coefficient  for  hemisphere 
cylinder 

is  a  body  station  non-dimensionalized 
with  respect  to  hemisphere  nose  diameter 


Equilibrium  gas  properties  at  the  stagnation  point  are  derived  from 
equilibrium  rormal  shock  tables  and  then  an  isentropic  expansion 
(i.e.  ,  body  streamline)  is  traced  on  a  Moliier  chart  to  the  calculated 
body  pressure.  Equilibrium  properties  are  then  read  off. 

Since  the  additional  variables  of  a  non-uniform  flow  field  with 
distance  from  the  nose  is  added  in  the  case  cf  a  blunt  body  a  specific 
flight  condition  has  been  chosen  as  an  example.  The  conditions  chosen 
are  V  =  25.000  ft/sec  and  h  =  250,000  ft.  The  surface  velocity  variation 
with  body  location,  X/D,  is  p.esented  in  Figure  3.19,  The  plasma  fre¬ 
quency  calculated  with  the  aid  ofReference  6  is  presented  inFigure  3.20 
along  with  the  plasma  frequency  and  collision  frequency  data  of  Reference 
7.  Tne  attenuation  as  a  function  ofX/Dis  presented  in  Figure  3.21.  Note 
that  for  the  example  chosen  an  antenna  would  have  to  be  located  about  100 
body  diameters  aft  of  the  nose  for  suitable  transmission  at  200  Me. 
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PLASMA 


Figure  3.1.  Bounded  plasma 
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everywhere  in  equilibrium. 

Same  as  EQ  but  with  boundary-layer  corrections  to 
temperature  and  density. 

Figure  3.  2.  Attenuation  on  a  blunt  body  versus  distance 
for  various  assumed  flow  fields. 
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ATTENUATION 


Figure  3.  3.  Variation  of  attenuation  with  plasma  frequency  ratio 
i  and  collision  frequency  ratio  for  a  bounded  plasma. 
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Figure  3.5.  Attenuation  as  a  function  of  pressure  and  temperature 
for  a  plasma  thickness  of  .2’. 
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Figure  3.  7.  Attenuation  as  a  function  of  pressure  and 
temperature  for  a  plasma  thickness  of  .  1'. 
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Cone  shock  wave  angle  versus  hyper 
similarity  parameter. 
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.15.  Surface  velocities  for  equilibrium  glide 
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Figure  3.  17. 


Plasma  and  collision  frequencies  versus  the 
equilibrium  glide  parameter. 
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Figure  3.  18.  Attenuation  versus  the  equilibrium  glide 
parameter  for  a  plasma  thickness  of 
.  1  ’  and  .  2' . 
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4.  USE  OF  JET  PENETRATIONS  AND  MIXING 
TO  CREATE  A  WINDOW 

4.  1  Basic  Principle 

One  of  the  most  attractive  methods  of  creating  a  "window"  in  a 
plasma  sheath  is  by  cooling  the  plasma  and  reducing  its  electron  con¬ 
centration  by  injecting  a  jet  of  electrophilic  gas  into  the  plasma.  The 
method  considered  herein  is  to  use  a  jet  from  a  flush  mounted  nozzle 
which  penetrates  into  the  flow  before  mixing  in  a  streamtube  carried 
downstream.  The  mixing  within  the  streamtube  serves  to  cool  the  plasma. 
Furthermore,  reduction  in  the  electron  concentration  by  attachment  to 
form  negative  ions  may  be  obtained  by  using  an  electrophilic  gas  as  the 
coolant. 

4.  2  Jet  Penetration 

In  order  to  cool  a  supersonic  flow  by  means  of  a  jet  the  maximum 
lateral  penetration  must  be  achieved  with  as  little  disturbance  of  the  main 
flow  as  possible.  The  specific  method  investigated  here  is  to  expandthe 
injected  fluid  to  a  high  Mach  number,  and  to  a  pressure  which  balances  the 
pressure  of  the  free  stream.  This  "correctly"  expanded  jet  is  directed 
across  the  main  stream  so  that  the  momentum  carries  the  jet  fluid  well 
out  from  the  orifice  before  being  swept  downstream.  The  fluid  spreads 
into  the  stream  by  turbulent  mixing.  The  behavior  of  the  flow,  the  mixing 
and  the  penetration  distances  will  depend  on  injection  Mach  number,  angle, 
fluid  density,  and  free  stream  properties. 

The  situation  is  shown  in  Figure  4.  1.  A  circular  jet  is  represent¬ 
ed  by  an  expanding,  curving  tube  in  the  free  stream.  The  tube  will  have 
more  and  more  free  stream  air  in  it  as  it  moves  outward;  this  air  is  en¬ 
trained  by  turbulent  mixing.  Near  the  orifice,  the  upstream  or  "leading" 
edge  will  act  like  a  bluff  obstacle  in  the  main  flow.  There  will  be  a  normal 
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shock  in  the  flow  before  it  spreads  around  the  jet  tube  and  closes  behind 
it.  In  reality,  part  of  the  jet  fluid  will  be  torn  away  from  the  tube  and 


_ _ i  j _ _ _ 

cl  x  i  ttru  uuwtiDt  x  cam, 
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total  mass  of  the  jet.  The  impact  of  the  free  stream  on  the  jet  tube  will 
cause  it  to  start  curving  downstream.  As  the  jet  progresses,  it  picks 
up  mass  by  entraining  free  stream  air,  while  maintaining  a  constant  mo¬ 
mentum.  Since  it  will  be  curving  there  will  be  a  centrifugal  force  which 
appears  in  the  force  equations.  As  the  jet  tube  turns  it  presents  a  dif¬ 
ferent  angle  to  the  free  stream,  thus  changing  the  impact  which  the  free 
stream  exerts  on  it. 


Two  torce  equations,  one  along  the  stream  tube  and  one  normal  to 
it  are  defined. 


a) 

b) 


aF  =  (  Am  )  V. 


aF  = 
n 


J 


{  Am  )  V.  tr 
J 


o  o 

m  c  AS 


(4.1) 


Am 


O 

m 


V. 

J 


AF 


AF 


n 


=  mass  of  jet  fluid  in  disc  of  diameter  W  ar.i 
thickness  S  along  the  jet. 

=  mass  flow  of  jet. 

=  average  velocity  of  jet  along  its  center  line 
at  station  S. 

=  approximate  path  length  from  orifice  along 
centerline  of  jet. 

=  arc  tan  (y/x). 

=  axial  force  acting  on  Am. 

=  normal  force  on  Am  {+  outward). 


Equations  (4.  1)  are  decoupled  by  the  following  assumptions.  First 
of  all,  the  jet  will  be  assumed  to  remain  intact  near  the  orifice  (i.  e.  the 
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jet  fluid  torn  away  will  be  neglected).  Second,  after  the  jet  has  made 
a  substantial  turn,  there  will  be  a  component  of  free  stream  velocity 
along  the  jet;  i  *  this  event  it  is  k*'Own  that  the  law  c:  entrainment  must 
be  modified^  {see  Pai,  1954  and  Townsend,  1956  ).  The  form  of 
modification  is  not  known  in  the  present  example,  and  will  be  neglected. 
Third,  the  impact  of  the  free  stream  takes  place  througl  a  normal 
shock,  so  that  the  pressure  surrounding  the  jet  near  the  orifice  will 
havv.  a  value  well  above  free  stream  static.  As  the  jet  turns,  the  shock 
becomes  oblique  and  the  pressure  rise  becomes  less.  This  decrement 
will  be  ignored  in  the  first  rough  force  balance.  The  effect  of  these  as¬ 
sumptions  is  to  fix  the  growth  of  mass  flow  in  the  jet  at  the  same  values 
as  would  occur  in  still  air.  Essentially  the  axial  force  balance  is  thus 
satisfied  by  hypothesis.  For  the  situation  ■where  the  stream  and  jet  are 

(9) 

different  fluids,  the  growth  rate  is  given  by  Ricou  and  Spalding  (1961)  . 

Their  formula  is: 


(4.  2) 


=  mass  flow  in  jet  at  station  S. 

2 

=  mass  flow  of  iet  at  orifice  =  (u/4  )  d  “  p.  V. 

'  'jo  jo 

=  density  of  air  stream  in  contact  writh  jet,  here 
taken  equal  to  ram  density  behind  a  normal 
shock.  It  is  assumed  that  ram  pressure  Pj^ 
of  the  c.ir  equals  the  jet  exit  pressure  (R)o. 

-  diameter  of  the  iet  orifice. 


Taking  the  momei  .am  in  the  axial  direction  as  constant,  thus  neglecting 
the  contributions  which  variations  :n  pressure  might  introduce,  and 
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neglecting  the  effects  ruled  out  by  the  second  a  isumption  above,  the 
average  velocity  of  the  jet  is  found  to  be. 


(4.3) 


The  width  of  the  circular  jet  W\  is  found  by  using  equations  (4.  2)  and 
(4.  3)  to  substitute  into  the  following  equation: 


m.  =  it/ 4  W.2  p  .  V.  (4.4) 

J  J  J  J 


Witn  p  .  taken  to  be  equal  to  since  the  density  of  the  jet  is 

J 

composed  largely  of  ram  air. 


Solution  for  W\  gives; 


W.  = 
J 


.  32  S 


(4.5) 


The  next  step  is  to  set  up  the  force  term  for  the  left  side  of  (4.1b); 
the  right  side  is  simply  the  centrifugal  force  of  acceleration.  The  force 
term  will  have  two  parts:  (1)  the  force  owing  to  the  accretion  of  free 
stream  air  to  the  jet;  and  (2)  the  force  owing  to  the  equivalent  aerodynamic 
drag  of  the  jet.  (Both  forces  cause  a  to  decrease  and  therefore  are  nega¬ 
tive).  The  first  of  these  depends  on  the  increase  in  mass  flow  per  unit  dis¬ 
tance  along  the  jet  dm/dS.  This  mass  before  uniting  with  the  jet  had  a 
velocity  normal  to  the  jet  of  U  °  Sin  <r  which  gives  rise  to  a  change  in 
momentum  normal  to  the  jet  or  a  normal  force  p6r  unit  length  of  jet  of: 

O 

=  -  U  °  Sin  <r  (4.  6) 

dS 

The  normal  force  due  to  aerodynamic  drag  is  given  by; 


1 

7 


U 


dA  CD  Sin 


(4.7) 
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According  to  Penland^*^,  the  drag  coefficient  Cn  is  constant  for  super¬ 
sonic  Mach  numbers  and  is  independent  of  or.  With  dA  =  WdS  and  W 
given  by  equation  4.  5  and  dm/  dS  obtained  from  the  Ricou,  Spalding  rela¬ 
tion  equation  4.  2,  the  force  equation  4.  1  can  be  written  as: 


(AF 

a5“ 


*1  <*) 


o  o 

-  —  m  cr 


(4.8) 


No'V  v/ith  substitution  of  equations  4.  6,  4.  7,  and  4.  2,  and  use  of  the 
relation 


•  dS  d<r  ..  do¬ 
or  =  ___  -  -Tf-r  "  V. 


■HT  H5 


j  HS' 


with  given  by  equation  4.3,  equation  4.  8  above  becomes: 


Sin  cr  + 


=  -  !  3 


2 

i*  VpjoAvjoApW 

•■(>)(£) 


D  >  S  *  Sin 2  0- 


1/2  1  dtr 


|H5“  *  d 


(4.  9) 


The  constant  in  the  second  bracket  of  equation  4.  9  will  be  treated  as  a 
scaling  parameter  in  the  following  way: 


1/2  p.  V  2V  p«  y/2 

ft  =  3.1  |  4-  1 

\Jzp.  U°2 


5ra  y 


=  3 


■<(>f 


1/2 
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where  the  dynamic  press- 
Newtonian  formula: 


at  the  jet  exit  is  given  by  the  .modified 


Pjo  '  PRA  =  !-9>  Si"  % 


By  defining. 


S  --  (S/d)  (  i  1 


a  -  2.94 


M.  sin  a 
, —  jo  o 

■ft  L - 7TTI72 

(p  RA/  p  j 


6  2- 


(fc) 


equation  ft  8  takes  the  form 


dtr  .  .  r  .  2 

—  r  —  sm  cr  —  kS  •  sin  <r 


(ft  10) 


The  equation  has  one  shape  parameter  k,  which  determines  the  precise 

form  of  the  function  cr  {§};  k  depends  only  on  free  stream  Mach  number. 

The  parameter  Q  contains  the  effect  of  jet  exit  Mach  number  and  angle  of 

injection.  Equation  4.10  cannGt  be  solved  analytically  but  can  be  easily 

(Hi 

programmed  on  a  Recomp  III  computer  as  described  by  Bartlett  .  The 
data  required  for  this  are  initial  conditions  and  values  of  k  and  a.  Once 
tr  is  obtained  as  a  function  of  S,  simple  quadrature  will  yield  the  co¬ 
ordinates  of  the  jet  centerline 


y/d  -  J"  sin  tr  (S,k)  d  (S/il)  -  a J' Bin  (r  (Sfk)  dS 
x/d  =  if  cos  <r  (S,  k)  d  S 


(4.  1  ila) 


{4.  1  lib) 


MC  6‘3-78-Rl 


For  the  computation  of  k  and  a  .  the  following  conditions  were 


chosen. 


V°  =  1.4 

if  =  14000  ft/  sec 
T*5  =  4 )  00°  K 
=  3.3 

PRA/P°=  4-U 


The  injected  fluid  has  a  Vj  =  1. 67  as  would  be  the  case  for  monatomic 
gas  in  the  jet; 


a  -  1.8M,  sin  cr 
jo  o 

k  =  .60 

The  solution  to  the  above  sets  of  conditions  are  shown  plotted  in 
Figures  4.  2,  4.  3,  4  4. 


4.  3  One -dimensional  Flow  with  Coolant  Injection 


4,  3,  1  Basic  Equations  lor  Perfect  Gases 

In  order  to  determine  properties  in  a  flow  field  after 
uniform  mixing  has  occured,  a  one-dimensional  stream  tube  analysis 
was  performed.  Schematically  the  problem  is  defined  in  the  sketch 
on  the  following  page: 
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Flow  at  (1)  is  known;  flow  properties  at  (2)  are  to  be  determined.  A 
coolant  or  electrophilic  gas  is  injected  into  the  stream  tube  between 
locations  (1)  and  (2)  with  known  properties.  Complete  mixing  takes 
place  in  the  region  between  (1)  and  (2),  and  the  mixing  of  the  initial 
flow  and  the  coolant  is  assumed  to  take  place  at  constant  pressure. 

Application  of  the  conservation  equations  for  mass;  partial  mass, 
momentum,  energy  and  a  state  equation  along  this  stream  tube  will 
determine  properties  of  the  mixed  flow'  at  (2). 

The  following  conservation  laws  result. 
plUlAl+mc  =  p2U2A2 
mc/plL’lAl  =  C2^  1  C2 

PjU^+m.U  =  P2U>2A2  (4.12) 

P1  U1  A1  H1  +  mc  Hc  =  p2U2A2H2 

P1  *1  T1  =  p2T2R2  Pl,  2 


39 


MC  63-78-R1 


- y-  - 


is  defined  as  the  percent  by  weight  of  coolant  present  in  the  mixed 
flow; 


P2(UC2)  =  "2..r 

P2(C2)  =  P2coolant 


(4.  13) 


With  this  definition,  weighted  values  for  Cp^  and  R^,  the  partial  en¬ 
thalpy  of  air  h^  and  of  the  coolant  h^  can  be  defined  as  follows, 
a 

h2  =  h2a  Jl  -  c 

CP2  ■  CPair  (‘‘C2)+CPc  H  <414) 

R2  *  Rair  (-2)  +  Rc  («2) 

With  assumption  that  the  moment  .m  of  the  injected  coolant  or  electro¬ 
philic  gas  is  much  less  than  the  momentum  of  the  incoming 

free  stream  p  ^  .there  results  after  defining  m/p^  Uj  Aj  =  Cq  , 

the  following  values  for  the  conditions  at  (Z). 

VU1  =  / 1  +  CQ, 

KljA  1  =  +  CQ^  '  P)/P2  (4-15) 

=  cq,  /*  +  ca, 

If  the  initial  conditions,  U^,  T^,  ^  etc.  are  speciiied  then  for  various 

values  of  the  parameter  CQ  cne  can  solve  for  the  downstream  conditions 

using  the  relations  above.  The  downstream  temperature  is  found  from  the 
energy  equation  to  be; 
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remains  constant  and  equal  to  its  free  stream  value  of 


F  =  h  + 
oe  cc 


(4.  17) 


An  alternate  approach  to  the  solution  of  this  problem  can  be  outlined  if 
one  wishes  to  determine  what  value  of  Cq  ,  and  U^,  resulted  in  a 
specified  downstream  temperature  for*a  certain  flight  velocity.  This 
is  done  in  the  following  way. 

Substitution  of  U_,  h^  from  relations  4.  15,  and  4.  14  into  the  energy 

equation  results  in  an  equation  for  U^. 


(4.18) 


Defining  the  stream  tube  mixing  problem  in  this  way  allows  us  to  select  the 
final  temperature  of  the  mixture  T^,  then  with  Hc  a  known  coolant  storage 
condition  and  h0  =  c  T_,  h,  =  c  .  T,,  U.  will  be  solely  a  function 

u  P  C-  L*  P  •  1 


:  T_,  h0  =  c 

p  L  2  p  . 

rc  a  *air 


of  Cn  or  the  ratio  of  coolant  to  initial  mass  flow  necessary  to  achieve  the 
U1 

desired  final  temperature. 


What  the  initial  temperature  must  have  been  in  order  to  reach  the  final 
temperature  with  the  proper  ratio  of  coolant  injected  into  the  flow 
can  be  found  from  equation  (4.  17)  to  be: 
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4.  3.  2  Special  Method  for  Real  Air  in  Thermodynamic  Equilibrium 


un - _ _ r _ 4. 

»»  UC  XI  v  1 1C  p?  1  xcv.  V  gc 


**««•»  «  «»  i  n  Z'  *i  *•  a!  <1.*^  V~»  1*^  f  V»  O 

U  O  D  UltipilV/ll  IP  V»»OVU*  vtvv*  vt«vr 


thermodynamic  properties  of  air  are  obtained  from  a  Mollier  chart 
(assuming  thermodynamic  equilibrium)  it  is  most  convenient  to  use 
the  method  in  which  the  downstream  temperature  is  specified.  It  is 
also  necessary  in  this  case  to  specify  the  local  pressure,  Pj  ^ 
flow. 


The  first  step  in  the  calculation  is  to  choose  a  value  of  the  partial 
pressure  of  air  at  station  2,  P^  ,  The  partial  pressure  of  the  coolant  is 
then 


P  =  P  _  p 
2c  ‘1,2  2a 


The  partial  densities  at  station  2  are 


Cdi 


2a 


R.  ZT_ 
A  2 


'2c 


P2c 

RcTZ 


assuming  Z  =  1  for  the  coolant.  The  mixing  ratio  is 


^2c  +  *^2a 

from  which  the  injection  coefficient  C~  can  be  computed  as 

U1 

CQt  =  C2/*-c2 

Since  the  partial  pressure  and  temperature  of  the  air  at  station 
2  have  been  specified,  its  enthalpy  is  also  determined.  Equation  4.  18 
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can  then  be  used  to  calculate  the  corresponding  value  of  U^,  the  local 
velocity  before  injection. 


This  straight  forward  method  yields  results  of  the  type  described 


in  the  next  section,  i.  e.  we  obtain  directly  curves  of  versus  U 


Q 


with  fixed. 


1 


r 


4.  3.  3  Selection  of  a  Coolant 

Since  data  is  yet  unavailable  to  compare  the  electrophilic 
effectiveness  of  various  gases  the  cooling  effectiveness  only  will  be 
considered  herein.  It  should  be  kept  in  mind  however  that  electrophilic 
additives  could  enhance  the  plasma  suppression  by  the  jet. 

The  effectiveness  of  a  coolant  may  be  judged  by  evaluating  the 
heat  it  can  absorb  as  it  is  raised  from  its  storage  temperature  to  the 
desired  temperature  of  the  coolant-air  stream  mixture  which  is  on  the 
order  of  3000°K  (5400°R).  As  may  be  recalled  from  Figures  3.4  through 
3.  7,  3000° K  would  lower  the  attenuation  of  a  200  Me  signal  through  a  0.  2 
ft  plasma  slab  to  about  3  db. 

Figure  4.  5  taken  from  Reference  12  shows  the  heat  capacity  per 
pound  of  hydrogen,  helium,  water  and  propane.  Hydrogen.the  most  effec¬ 
tive  coolant  could  only  be  used  where  it  could  be  shown  that  combustion 
would  not  occur,  or  if  it  does,  would  not  be  harmful.  Water  was  included 
for  comparison  as  NASA  has  been  engaged  in  an  intensive  effort  in  plasma 
suppression  by  water  injection.'  However,  the  prime  interest  in  the  studies 
at  hand  is  in  a  light  gas  system.  Helium  and  propane  are  quite  comparable 
in  the  5000° R  „*ange  of  interest. 

Specific  studies  in  section  4.  3.4  will  be  based  upon  helium  as  the 
coolant. 
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4.3.4  Results  of  Coolant  Injection 

The  equations  developed  in  section  4.  3.  2  allow  us  to 
calculate  an  injection  coefficient: 

CQi  =  “/P!  U1  A1 

to  obtain  any  desired  temperature  at  station  2  by  injecting  a  coolant  into 
a  given  stream  at  station  1.  Since  the  desired  temperatures  cover  a  wide 
range  and  the  stream  at  station  1  can  correspond  to  many  different  antenna 
locations  on  vehicles  travelling  at  many  speeds,  a  concise  presentation  of 
the  results  of  our  calculations  requires  some  unifying  scheme. 

We  have  chosen  to  consider  the  problem  in  the  following  terms; 

a)  the  basic  independent  parameter  is  the  vehicle 
flight  speed;  together  with  nominal  atmospheric 
temperature  of  41 9° R,  this  defines  a  stagnation 
enthalpy  H^. 

b)  antenna  location  on  the  vehicle  is  defined  in  terms 

of  the  local  velocity  of  the  stream  external  to 

it.  For  antennas  located  on  slender  vehicles  or  on 

the  shoulder  cf  blunt -nosed  vehicle  U,  0.  6  U 

1  o c 

well  beyond  the  sonic  point. 

c)  the  coolant  temperature  is  fixed  at  540*R. 

d)  the  coolant  is  assumed  te  be  helium,  chosen  as  a 
typical  example  of  a  coolant  with  a  large  capacity 
to  absorb  heat. 

The  energy  equations  4. 16  and  4. 18  show  that  it  is  possible  with 
this  scheme  to  either 

1)  cho'-se  U.  and  obtain  a  relation  between  C-.  and  T-. 

1  Q,  2 

r 

2}  choose  a  desire*.  and  obtain  a  relation  *  i'veen 
CQ  and  W 
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jooth  methods  of  presentation  are  used  for  the  results  shown  in 

Figures  4.6,  4.7,  4.8,  4.9.  Rather  than  use  Cq  ,  however,  which  is 

U1 

normalized  with  respect  to  a  streamtube  area  at  the  injection  station  we 
have  used 


m 

c 


P1  °!A 


2 


which  is  based  on  the  streamtube  properties  before  injection  and  the  area 
actually  cooled  at  the  antenna  stat  .on.  This  parameter  can  be  related 
more  directly  to  the  antenna  requirements. 

Most  of  the  calculations  so  far  have  been  performed  for  the  sake 
of  simplicity  v/ith  the  assumption  that  air  behaves  as  a  perfect  gas.  One 
set  has  been  performed  using  equilibrium  air  properties  to  illustrate  the 
differences  to  be  expected  in  general.  The  coolant  in  all  cases  has  been 
helium. 

Figures  4,  6,  4.  7,  4.  8  illustrate  how  the  injection  requirement 
varies,  with  the  desired  end  temperature.  Three  flight  speeds  are  repre¬ 
sented  and  curves  are  drawn  for  constant  values  of  the  local  velocity  Uj. 

Note  that  each  curve  exhibits  a  peak  at  a  value  of  T^»  which  lies  in  the 
range  of  interest  for  plasma  suppres  sion.  Thus  for  a  given  flight  speed  and 
local  velocity  the  injection  squirem  *.nt  is  not  very  sensitive  to  the  speci¬ 
fication  of 

Figure  4.  9  is  a  cross  plot  of  ihese  data  for  T^  =  2400°K  =  4325*R. 

It  illustrates  how  the  coolant  requirement  falls  off  as  the  antenna  station  is 
moved  downstream  to  higher  values  of  U^.  The  curve  for  a  flight  speed  of 
18,  000  ft/sec  is  reproduced  in  figure  4. 10  where  it  is  compared  with 
equilibrium- air  results  for  three  values  of  the  local  pressure.  The  injection 
requirement  predictions  are  considerably  lower  for  equilibrium  air  than  for 
air  behaving  as  a  perfect  gas. 

The  shaded  zones  in  the  lower  right-hand  corners  of  some  of  tfhe 
figures  correspond  to  heating  of  the  stream  rather  than  cooling.  The  *sxr 
pianation  for  this  unexpected  behavior  lies  in  the  fact  that  coolant  injection 
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decreases  the  velocity  of  the  airstream  as  well  as  reducing  the  total 
enthalpy  oi‘  the  mixture.  The  decrease  in  velocity  raises  the  ratio  of 
static  temperature  to  total  temperature  and  under  some  conditions  this 
process  outweighs  the  total -temperature  decrease.  The  practical  con¬ 
sequence  is  that  injection  rates  must  exceed  a  certain  critical  value 
before  any  cooling  is  accomplished. 


4.4  Effect  on  Vehicle  Aerodynamics 

The  use  of  a.  side  jet  to  create  a  "window"  will  create  a  moment, 
upsetting  the  vehicle  trim.  In  order  to  get  an  order  of  magnitude  esti¬ 
mate  for  trim  angle  with  the  side  jet,  a  Trailblazer-type  vehicle  with  the 
characteristics  shown  in  Figure  4.11  will  be  considered.  The  flight 
conditions  will  be  assumed  constant  as  follows: 


U  =  18,  000  ft/  sec 

oc  ' 

T  =  233°K 

oc 

VU«  =  0.67 

pl/p«c  =  0.52 


The  coolant  will  be  aken  as  helium  and  has  a  m 
as  shown  in  Figure  4 JO  for  an  equilibrium  gas. 

For  a  jet  penetration  into  the  shock  layer,  the  axial  force  normal 
to  the  surface  of  the  missile,  N^,  has  been  assumed  to  be  twice  the  exit 
jet  momentum  or. 


/?1  Uj  A2  (C0  )  of  0.13 


N.  =  2  mU. 

J  J 


(4.19) 


but 


m 


:Q?  Pi*iA. 
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so 


Nj  =  0.26  p  j  Uj  A2M^.  a^ 


(4.  20) 


*6 


The  area  A ^  is  the  area  that  the  jet  cools  and  will  be  taken  as 

where  6  is  the  plasma  thickness.  From  Figure  4. 11  we  find  that 
2 

A 2=  0,  001 96f  .  Also  from  Figure  4. 11  note  that  the  trim  angle  of 
attack  is  given  by: 


a 


Trim 


127.6  N. 
_ _  J 

“  T 


(4.  21) 


CN  V 
a 


Substituting  equation  (4.  20)  and  the  above  value  of  into  (4.  21) 
we  find 


Trim 


.  °*13  Pi/Pce  V^Vj 


(4.  22) 


'N  U 
a  c 


It  now  remains  to  find  the  jet  Mach  number,  M..  In  order  to 

J  _ 

cool  the  pUvpmit  layer  of  thickness,  6,  it  will  be  necessary  that  the  , 
jet  penetrate  a  distance  of  6/2.  For  a  typical  jet  the  maximum  jet 
penetration  was  given  in  Reference  13  as 

y  /d.  =  2.  2  Mt  (4.  23) 

'or/  j  J  .. 

e  said  that  we  war 

found  from  the  mass  flow  of  coolant  required. 


We  have  said  that  we  want  =  6/2.  The  jet  exit  diameter  maybe 


m  =  *p  ,,  M.  a.  A  -  =  PT  MT  A. 


J  J  J 


J  J 


N 


irjT7 


ft 

m 


But  from  section  4.  2, 


Pj  =  1. 8q  =  0.9  P1U1< 


therefore 


m.  =  0.  9p,  U.2M.A> 
3  1  3  J 


M 


vi 


R*  T* 
J  J 
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As  stated  before  m  =  C_  p,  U,  A. 

U0  1  1  t 


:Q2  p  1  U1  A2  =  °*  9  Pi  U1  Mj  Aj  ^  Vj  /  RJ  T J 


0.9  U.Mt 

a2/aj  =  .-!_! 


j  *  j 


(4.  24) 


where 


AO 

2  4.0 


d. 

A.  =  7T — — 

J  4 


therefore 


S2  _  0-9U,Mj  |  Yj 


\~*P7 


From  (  4.23) 


(^7) 


\ 

21.  5  Cn 


(4.  25) 


or  for  the  values  in  the  example 


Mj  ^  2.16 
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The  trim  angle  of  attack  may  now  be  found  by  substituting  in  equation(4.  22) 


a 


Trim 


(0.13)  (0.52)  (0.67)  (2. 16)  (3330) 
(.03)  (18,000) 


a 


Trim 


0.  6" 


Such  a  value  for  the  trim  angle  of  attack  is  not  considered  a  problem 
of  much  concern.  In  fact  should  a  symmetrical  installation  be  used  the  trim 
angle  would  be  zero. 
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Figure  4.4.  Jet  centerline  for  M*  =  3  —  4  and  various  jet 
mach  numbers. 
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Figure  4.  7 


Tt~»*  *K 

Coola.nt  requirement  as  a  function  of  temperature 
for  a  free  stream  velocity  of  21,000  ft/sec. 
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5.  USE  OF  A  FLOW  DIVERTER  AND  BASE  COOLING 
TO  CREATE  A  WINDOW 


5.  J  Basic  Principles 

Re-entry  vehicles  are  enveloped  by  a  thin  region  of  ionized  gases 
which  hinder  communication  from  the  vehicle  to  monitoring  stations  on  the 
ground.  A  method  of  alleviating  this  communication  problem  would  be  to 
somehow  eliminate  this  region  of  ionized  flow  locally  over  an  antenna 
surface,  thus  creating  a  window  through  which  transmission  is  possible. 
Locating  a  flow  diverter  just  upstream  of  an  antenna  location  is  an  attempt 
t.o  create  such  a  window  by  diverting  the  flow  of  ionized  gases  away  from 
the  antenna  surface. 

At  the  base  of  the  flow  diverter  a  "dead  air"  region  will  exist  which 
will  contain  temperatures  equivalent  to  those  in  a  stagnation  point.  There¬ 
fore  a  coolant  is  injected  into  this  region  to  keep  temperatures  here  at  a 
tolerable  level  for  communication. 

The  simplest  flow  di-  erter  to  consider  is  one  which  extends  outfrom 
the  vehicle  surface  perpendicular  to  the  stream.  All  the  oncoming  flow 
would  be  diverted  laterally  and  a  window  would  be  created  at  its  base  (Figure 
5.1  A),  Unfortunately  ,  such  a  protruberance  would  be  bard  to  protect  against 
aerodynamic  heating  and  would  include  significant  aerodynamic  forces  on  the 
vehicle.  Generally,  these  considerations  would  eliminate  the  practical  use  of 
this  type  of  two  dimensional  diverter  with  the  exception  of  lifting  re-entry  ve¬ 
hicles  with  a  "natural"  flow  diverter  such  as  a  dorsal  fin. 

A  more  versatile  diverter  must  have  a  swept-hack  leading  edge  and  be 
shaped  more  like  a  cone  or  pyramid  (Figure  5.  IB).  Now,  not  all  the  flow  is 
diverted  around  the  side  and  some  flows  over  the  top  for  the  antenna  to  look 
through.  This  kind  of  flow  diverter  will  be  useful  only  if  the  flow  over  the 
top  is  much  more  transparent  than  the  original  flow.  Of  course,  the  use  of  a 
diverter  is  only  feasible  in  the  first  place  if  the  region  of  high  attenuation  is 
more  or  less  concentrated  in  a  finite  layer  near  the  surface. 
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5.2  Diverter  Effectiveness 


For  a  communication  gain  to  be  realized  with  a  flow  diverter  the 
plasma  sheath  must  be  either 

a.  completely  diverted  from  flowing  over  the  antenna 
at  the  base  of  the  flow  diverter. 

b.  geometrically  thinned  such  that  the  plasma  slab 
thickness  at  the  antenna  location  is  significantly 
less  than  the  initial  plasma  thickness. 

Since  two  dimensional  flow  diverters  are  only  practical  under 
special  circumstances,  due  to  drag  and  heating  penalties,  we  will  concen¬ 
trate  on  the  three  dimensional  flow  diverter  which  thins  the  plasma  sheath. 
Figure  5.  2  shows  the  decrease  in  attenuation  with  decreasing  plasma 
thickness  for  an  antenna  located  on  the  aft  portion  of  the  conical  surface  of 
a  Trailblazer  II  vehicle  re-entering  the  atmosphere  at  18,  000  fps.  A.  homo¬ 
geneous  finite  thickness  slab  model  of  the  plasma  has  been  used  for  these 
calculations,  and  the  local  flow  has  been  assumed  to  be  in  equilibrium. 

As  a  means  of  geometrically  thinning  the  plasma  sheath,  2  types 
of  flow  diverters  have  been  studied 

a.  conical  flow  diverters 

b.  conical  flow  diverters  of  elliptic  cross-section 

The  elliptic  cones  were  analyzed  since  it  was  felt  that  the  cross-flow 
associated  with  this  shape  would  contribute  to  the  thinning  of  the  plasma 
sheath  by  diverting  part  of  the  hot  ionized  air  to  the  sides  of  an  antenna. 

In  the  limit  of  very  large  elliptic  ratios  the  flow  over  an  elliptic  cone  flow 
diverter  is  seen  to  act  like  a  two  dimensional  flow  diverter. 

5.  2. 1  Conical  Flow  Diverter 

Figure  5.  3  is  a  schematic  illustration  of  the  flow  past  a 
diverter  which  consists  of  half  a  circular  cone.  The  flow  approaching 
the  diverter  is  first  compressed  in  flowing  over  the  conical  surface  and 
then  re-expanded  as  it  leaves  the  base.  The  flow  over  the  antenna  is  thus 


62 


MC  63-78-R1 


i 


analogous  to  that  .n  the  near  wake  of  a  full  cone.  If  we  assume  that  the 
pressure  at  the  antenna  station  is  roughly  equal  to  that  in  the  undisturbed 
stream,  we  can  calculate  the  thinning  of  the  high-attenuation  layer  in  the 
plane  of  symmetry  as  a  purely  geometrical  effect.  If 

6j  is  the  original  thickness  the  high  attenuator  layer 
6^  is  its  thickness  over  the  antenna 

R  is  the  radius  of  the  base  of  the  diverter 

conservation  of  streamtube  area  requires  that 
52  =  4^777  -  R 

This  relationship  is  illustrated  in  Figure  5. 4  by  a  plot  of  6~/5.  versus 

R/6r 

The  resulting  gain  in  signal  strength  is  illustrated  in  Figure  5.  5 
which  combines  the  results  in  Figure  5.2  and  5.4  to  show  the  decreased 
attenuation  versus  conical  ixow  diverter  radius. 

Bearing  in  mind  the  3  db  criteria  assumed  for  adequate  transmission 
the  reduction  in  attenuation  is  obviously  inadequate.  Figure  5.  6  is  merely 
an  extension  of  Figu  re  5.  5  to  determine  just  what  flow  diverter  heights  are 
required  to  meet  thi  5  criteria.  Note  that  the  flow  diverter  heights  required 

J 

‘  are  prohibitive. 


% 


iT7 


5.  2.  2  Elliptic  Cone  Flow  Diverter 

An  elliptic  cone,  as  shown  in  Figure  5.7  experiences  the 
same  geometrical  thinning  of  the  plasma  sheath  as  the  circular  cone.  As 
before  conservation  of  streamtube  area  requires  that 


b,  a,  major  and  minor  axes  of  the  ellipse. 


■ 
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This  relationship  is  exactly  the  same  as  that  shown  plotted  in  Figure  5.4  if 


In  order  to  make  a  comparison  of  the  relative  thinning  effects  of  the 
elliptic  and  circular  base  cone,  we  define  an  equivalent  circular  cone 
as  one  whose  base  area  is  equal  to  that  of  the  elliptic  cone.  This  de¬ 
termines  an  equivalent  radius  for  the  circular  cone  to  be  the  following 

^equivalent  ~  >ja 

To  determine  the  relative  merit  of  the  circular  and  elliptic  cone  in 
thinning  the  high  attenuation  layer  for  the  conditions  shown  below. 


we  have  to  refer  to  Figure  5.  5  and  note  the  attenuation  for  various 

a4>K 

altitudes  for  R  .  =  1.414  and  R  = — - —  =1.5. 

equivalent  2 

The  elliptic  cone  decreases  the  attenuation  of  a  signal  at  ail  altitudes 
by  approximately  .5  db  more  than  the  circular  cone.  Hardly  significant 
in  comparison  to  the  absolute  value  of  the  attenuated  signal. 
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Howeyfi?  by  mounting  an  elliptic  flow  diverts?  with  its  major 
axis  perpendicular  to  the  flow  we.  can  take  advantage  of  the  cross - 
flow  whose  direction  is  from  the  major  to  the  minor  axis  to  divert 
part  of  the  plasma  layer  to  the  sides  and  thereby  further  thin  the 
plasma  layer  over  an  antenna. 

Data  taken  from  References  14,  15,  16  has  been  used  to  calcu¬ 
late  the  shock  shape  about  two  e'liptic  cones  with  }xis  ratios  of  1.  394 
and  1. 788.  The  Mach  number  of  the  flow  over  these  cones  was  taken 
to  be  M  -  3.  088  and  M  =  5.  9409  for  each  of  the  two  elliptic  cones. 

Figure  5.  8  shows  the  elliptic  cones  used  in  the  analysis  and 
their  shock  envelopes.  The  equivalent  circular  cone  as  defined  pre¬ 
viously  above  is  that  circular  cone  whose  cross-sectional  base  area 
is  equal  to  that  of  the  elliptic  cone.  The  antenna  looks  out  through  a 
window  vhose  boundary  is  shown  by  the  dotted  lines. 

As  a  means  of  comparing  the  effects  of  how  thick  the  shock 
layers  are  for  both  the  equivalent  and  elliptic  cone  flow  diverters.,  the 
ratio 

Ab  ...  /aR 
ellipse/  cone 

has  been  plotted  in  Figure  5.  9  versus  elliptic  axis  ratio  b/a  for  constant 
Mach  numbers.  The  mass  flow  contained  in  the  shock  layers  for  both 
the  conical  and  elliptic  flow  diverters  are  approximately  equal,  so  that 
the  cross-flow  effect  in  diverting  the  flow  to  the  sides  of  an  antenna  can 
be  determined  from  this  figure.  It  is  seen  that  for  a  Mach  number  3.  0 
which  corresponds  to  the  local  flow  velocity  at  possible  antenna  locations 
on  the  Trailbiazer  G  vehicle,  the  2-dimensional  effect  of  diverting  flow  to 
the  sides  with  elliptic  cone  flow  diverters  is  sign  Jicant.  For  a  cone  of 
elliptic  axis  ratio  of  1. 5  the  plasma  thickness  is  3/4  that  of  the  flow  over 
the  equivalent  circular  cone.  Referring  to  Figures  5.2  and  5.4  with  "a* 
assumed  equal  to  1. 0  and  the  initial  plasma  thickness  =  2*  the  decrease 
in  attenuation  due  to  the  cross-flow  associated  with  the  elliptic  cone  is 
approximately  3  db  greater  than  that  derived  by  use  of  the  equivalent  cir¬ 
cular  cone.  This  is  again  a  relatively  low  gain  in  comparison  to  the  ' 
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absolute  magnitude  of  the  signal  loss.  However  it  does  point  out  the 
effect  of  cross -flow  on  thinning  the  plasma  sheath  in  comparison  to 
solely  geometrical  thinning. 

% 

5.  3  Cooling  the  Near  Wake 

Even  for  the  above  gains  to  be  realized  it  is  of  course  necessary 
to  maintain  low  electron  concentrations  in  the  separated  flow  behind  the 
base  of  the  diverter.  This  is  accomplished  by  injecting  a  coolant,  such 
as  helium,  into  the  base  flow  so  as  to  avoid  recirculation  of  the  hot  external 
flow. 

The  flow  diverter  scheme  as  a  method  of  creating  an  antenna 

(17) 

window  is  sketched  in  Figure  5. 10.  As  Chapman'  points  out,  if  gas 
is  injected  at  least  as  fast  as  the  mixing  layer  scavenges  it  from  the 
dead-air  region,  no  flow'  is  recirculated  at  the  base,  and  the  temperature 
in  the  dead-air  region  is  equal  to  that  of  the  injected  gas.  It  does  not 
matter  whether  the  mixing  layer  is  laminar,  as  considered  in  some  detail 
by  Chapman,  or  turbulent;  however,  the  rate  of  gas  injection  required  to 
avoid  recirculation  is  much  larger  in  the  latter  case.  The  theory  pre¬ 
sented  here  will  concern  itself  with  turbulent  flows,  which  will  occur  in 
most  cases  of  practical  interest.  In  particular  we  consider  rates  of  injec¬ 
tion  less  than  those  required  to  avoid  recirculation  and  calculate  the  varia¬ 
tion  with  injection  rate  of  the  temperature  in  the  dead-air  region. 

Since  very  little  information  is  available  about  the  structure  of 
turbulent  layers  in  high-speed  flows,  a  simple  analysis  seems  most 
appropriate.  We  shall  assume  that  the  profiles  of  total  enthalpy  and 
foreign  gas  concentration  are  linear  functions  of  the  velocity  profile. 


corresponding  to  =  L ,  =1.0 

II  =  H  +  (H  -  H  )  u* 
s  e  s' 

(5.1) 

C  =  C  ( 1  —  u*  ) 
s  '  ' 

(5.  Z ) 
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is  the  local  stagnation  enthalpy 

is  the  enthalpy  of  the  separated  region 

is  the  stagnation  enthalpy  of  the  free  stream 

is  the  local  foreign  gas  concentration 

is  the  foreign  gas  concentration  in  the  separated 

region 

is  the  local  velocity  divided  by  that  of  the  free 
stream. 


We  shall  also  assume  the  form  of  the  eddy  viscosity  law  sug~ 

(1 81 

gested  by  Ting  and  Libby  ,  which  relates  the  compressible  profiles 
to  the  incompressible  profile  through  the  Howarth  transformation: 

Y  =  /  p  *  dy  (5.3) 

At  reattachment  the  flow  in  the  mixing  layer  splits  as  sketched 
in  Figure  5. 10.  The  dividing  streamline  (  y=  0)  splits  the  flow  such 
that  the  mass  flux  flowing  below  it  is  equal  to  the  total  rate  of  mass  en¬ 
trainment  from  the  dead-air  region^1 7^.  The  reattachment  streamline 
{v  =  y-  )  splits  the  flow  such  that  the  mass  flux  flowing  below  it  is  recir¬ 
culated. 

if  b  is  the  width  of  the  flow  (or  the  perimeter  in  the  axisymmetric 
case),  the  following  relation  holds: 

/c 

p  udy 

id. 


the  entrainment  of  foreign  gas:  bC 


f%  u, 

J  -oo 


ryi 

the  recirculation  of  foreign  gas:  b  f  p  uCdy 
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In  terms  of  normalised  variables,  the  balance  of  in-flow  and  out- flow 
of  foreign  gas  may  thus  be  written: 


C 

s 


p  *  u*  dy 


p*  u*  dy  +  Cs 


p  *  u *  C  *  dy  (5.4) 


Using  (5.  2)  and  (5.  3)  this  may  be  solved  for 


C 

s 


1 

ITT 


(6.  5) 


where 


(5.  6) 


For  the  case  of  zero  heat  transfer  to  the  solid  surface  bounding 
the  separated  region  (which  will  yield  a  "recovery  factor"  with  injection), 
a  similar  balance  may  be  written  for  the  flux  of  total  enthalpy: 


p  *  u  *  dy 


jr°  ,yi 

r  p  *  u*  dy  +  C  p 

yi  •/-» 


Using  (5.1),  this  reduces  to: 


J 

J  T  1 


i  -  refers  to  injected  coolant  conditions. 


u*Hdy  (5.7) 


(5.  8) 


For  perfect  gases  with  constant  specific  heats,  this  may  be  expressed 
in  terms  of  temperature  as: 


J 


j  + 


Pi 


(5.  9) 
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New  J  may  be  related  to  the  location  of  ine  re -attachment  stream¬ 
line  in  the  flow  by  (5.  6),  which  is  a  transformed  expression  to  be  evaluated 
from  the  solution  for  incompressible  flow.  The  same  solution  can  be  used 
to  evaluate,  the  ratio  of  the  injection  rate  to  the  total  entrainment  rate 
of  the  mixing  layer: 


(5.10) 


The  relation  between  4>  and  J  is  given  in  Figure  5. 11  with  the  assumption 
made  here  that  it  is  universal  to  all  turbulent  mixing  layers. 

(19) 

Hill  and  Nicholson'  '  have  measured  turbulent  entrainment  rates 
over  a  very  wide  range  of  conditions  and  have  developed  the  following 
formula  for  the  entrainment  coefficient  c  : 

q 


P  u 
r  s  s 

p  u 
r  e  e 


=  0. 049 


r  i 

R  T 

0.4 

’  T 

e  e 

e 

TTT — 
8  8 

^e 

m  •» 

1  0.67 


(5.11) 


In  this  expression  the  gas  constant  ratio,  R  /  R  ,  in  the  separated  region 

8  6 

may  be  evaluated  in  terms  of  the  molecular  weight  ratio,  m  /  m.,  and  the 

-C  X 

mixing  ratio  Cs: 


R  /R  =1.  +C 
s'  e  s 


(5.12) 


This  coefficient  c  has  been  normalized  with  respect  to  the  base  area  of 

q 

the  flow  diverter  (Reference  20)  such  that  proper  comparison  can  be  made 
with  the  coolant  injection  jet-penetration  coefficient.  It  was  found  that; 


o 

m. 

i 


pl  U1  a2 


6  c 

q 


To  determine  the  coolant  mass  flux  CQfor  varying  temperatures  iu  the 
mixing  layer  greater  than  the  injected  coolant  temperature, equation  5.  9 
is  solved  for  the  parameter  J  which  from  figure  5. 11  defines  the  ratio 
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c.f  injection  rate  to  the  total  entrainment  rate  of  the  mixing  layer. 

Solution  of  equation  5.  11  defines  the  total  entrainment  of  the  mixing 
layer.  Thus  curves  showing  the  effect  of  mass  injection  as  a  function 
of  separated  region  temperature  can  be  drawn. 

Figure  5. 12  shows  Cq  plotted  versus  local  flow  diverter  velocity 
or  body  station  for  various  separated  region  temperatures. 

The  temperature  in  the  mixing  layer  in  general  varies  mono- 
tonically  from  T  ,  the  temperature  in  the  separated  region,  to  T  , 
the  temperature  in  the  free  stream  outside  the  region  of  influence  of 
the  flow  diverter.  However  for  some  combinations  of  T  and  the  local 
flow  velocity  Ue>  at  the  flow  diverter  station,  the  mixing  layer  temperature 
varies  nonmonotonic  ally  from  its  two  boundary  conditions,  and  in  this  case 
a  peak  temperature  can  exist  which  is  higher  than  the  free  stream  tempera¬ 
ture  T  .  When  these  conditions,  shown  circled  in  Figure  5.  12,  occur, 
e 

there  is  obviously  no  gain  in  cooling  the  local  flow  by  means  of  a  flow 
diverter.  This  effectively  sets  a  limit  on  flow  diverter  location. 
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Figure  5.5.  Attenuation  versus  flow  diverter  height 
as  a  function  of  altitude. 
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Figure  5.7.  Schematic  of  elliptic  cone  flow 
diverter  performance. 
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6.  COMPARATIVE  EVALUATION  OF  THE  TWO  WINDOWS 


The  work  reported  in  the  previous  sections  has  been  concerned 
with  two  methods  of  modifying  and  cooling  the  flow  field  in  the  vicinity 
of  an  antenna  installation  on  a  re-entry  vehicle. 

In  one  method  nozzles  flush  with  the  vehicle  surface  are  used  to 
project  jets  of  coolant  into  the  flow  .  Ac  it  is  carried  downstream,  the 
coolant  cools  the  airstream  by  mixing  with  it.  In  the  other  method  a 
physical  flow  diverter  is  used  to  displace  the  highly  ionized  flow  laterally 
around  the  antenna  aperture,  A  region  of  separated  flow  is  formed  over 
this  aperture  behind  the  base  of  the  flow  diverter.  Coolant  is  injected  to 
maintain  low  temperature  both  in  this  region  and  in  the  shear  layer  which 
bounds  it. 

The  injection  of  a  coolant  with  a  flush  mounted  nozzle  appears  to 
be  by  far  the  most  attractive  and  practical  method  for  the  following  reasons: 

1.  The  jet  method  can  reduce  the  plasma  temperature  to 

levels  at  which  transmission  is  possible  (attenuation 
less  than  3  db)  with  modest  mass  flows  (see  Figure  4. 10}. 


The  flow  diverter  method  although  reducing  the  attenuation  will  not 
reduce  the  attenuation  to  the  3  db  level  as  was  demonstrated  in  Figures 
5.  5  and  5.  6. 

2.  The  jet  method  allows  the  cooling  cf  a  plasma  sheath 
located  away  from  the  body  surface  whereas  the  flow 
diverter  is  limited  to  relatively  thin  plasma  sheath 
near  the  body  surface. 

3.  The  jet  method  allows  all  equipment  to  be  located  in¬ 
ternally  while  the  flew  diverter  is  an  external  pro¬ 
tuberance  which  must  be  protected  against  aero¬ 
dynamic  heating  and  loads. 
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